jVteA-eJt-JWjSOV 


NASA Contractor Report 172307 


^ASA-CR ' 172307 

-^ 93400128 ^^ 


STRESS INTENSITY FACTORS IN A REINFORCED 
THICK-WALLED CYLINDER 


Renji Tang and F. Erdogan 


LEHIGH UNIVERSITY 
Bethlehem, Pennsylvania 18015 


Grant NGR 39-007-011 
February 1984 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 




STRESS INTENSITY FACTORS IN A 
REINFORCED THICK-WALLED CYLINDER* 


by 

Renji Tang and F. Erdogan 
Lehigh University, Bethlehem, PA 


Abstract 

In this paper an elastic thick-walled cylinder containing a radial 
crack is considered. It is assumed that the cylinder is reinforced by an 
elastic membrane on its inner surface. The model is intended to simulate 
pressure vessels with cladding. The formulation of the problem is reduced 
to a singular integral equation. Various special cases including that 
of a crack terminating at the cylinder-reinforcement interface are inves- 
tigated and numerical examples are given. Among the interesting results 
found one may mention the following: In the case of the crack touching 

the interface the crack surface displacement derivative is finite and 
consequently the stress state around the corresponding crack tip is 
bounded, and generally, for realistic values of the stiffness parameter, 
the effect of the reinforcement is not very significant. 


1. Introduction 


Because of its practical importance the problem of a thick-walled 
cylinder containing a radial crack has been studied rather extensively 
(see [1] for review and references). Considered as a plane problem 
the results obtained from such a study provide an upper bound for the 
stress intensity factors in pressure vessels, pipes, and disks containing 
an axial part-through crack. In all these previous studies it is assumed 
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that the cylinder is a homogeneous elastic solid. However, there is 
at least one group of technologically important problems in which the 
material may not be considered homogeneous. These are the problems 
Which may generally be classified under "multi-layered" cylinders. 

Variety of reinforced disks and cylinders and pressure vessels with 
cladding may be cited as some examples, in this paper we are primarily 
interested in thick -walled cylinders with cladding containing a radial 
Crack. This is a nonhomogeneous cylinder problem the distinguishing 
feature of which is that the thickness of the clad (that is, the inner 
cylinder) is very small compared to that of the main cylinder. The 
crack problem may be formulated and solved directly as a nonhomogeneous 
dylinder problem. However, there are at least two reasons for not 
following this procedure. First, the analysis of the nonhomogeneous 
problem is highly complicated and secondly if, as in the cladded pressure 
vessels, the thickness ratio of the two cylinders is very small , past 
experience i ndi cate that there would be severe convergence problems in 
the related numerical calculations. In this paper the cladding will 
instead be approximated by a membrane and the problem will be treated 
as a reinforced thick -walled cylinder. A similar problem for a half 
plane and an infinite strip is considered in [2] and [3]. 

Introducing a dislocation into Mitchell's solution the problem 

1 

is formulated in terms of a singular integral equation. The case of the 
crack terminating at the interface is separately studied and the stress 
intensity factors for various values of the stiffness constant character- 
izing the reinforcement are given. 


2. The Basic Formulation 

Except for the application of the boundary conditions, the formulation 
of the problem is quite similar to that given in [1]*. The stress state 
in the reinforced cylinder will again be expressed as the sum of two 


Pi 


In this paper certain expressions have been simplified through fur- 
ther manipulations and by evaluating some of the series in closed 
form. 



solutions containing the main features of the crack and the hollow 
cylinder, respectively. Omitting the details and referring to Figure 
1 these two solutions may be written as follows (see [1] and [4]-[6]): 
a. The cAXLck iohxtion 


a, (r,e) = - --fow f r- r cose-t-2t sin 2 e _ 2t 2 sin 2 e(r cose-t) 1 r , n „ 
lrr it(k+ 1) J 1 *r4t z -2rt cose (r 2 +t 2 -2rt cose) 2 


( 2 . 1 ) 


«. L f r. s i ne(2t cos9-r) 2t sine(r cose-t)(r-t cose 1 v .. 

°lre^ r ’^ ^+ry J L ?2+t2-2rt cose F n^-Irt cose) 2 L Jf(t)dt, 


( 2 . 2 ) 



2 u 

TT ( K+l ) 


r 2cose(r-t cose)+r cose-t 
L r 2 +t 2 -2rt cose 


2(r cose-t)(r-t cose) 2 
(r z +t 2 -2rt cose) 2 t_ 


]f(t)dt , 


(2.3) 


where y is the shear modulus, k = 3-4v for plane strain, k = (3-v)/(l+ v ) 
for plane stress, v the Poisson's ratio and the function f is defined by 

f(r) * [u Q (r,+0) - u Q (r,-0)] , (e<r<g) . (2.4) 

In order to facilitate the application of the boundary conditions, 
it is convenient to express the stress components on the boundary in 
terms of the following Fourier series: 

r 9 ~ 

a lre (a,0) = ¥T^+TT J ?, A n (t)stn ne » (2.5) 

e n "‘ 

o lrr (a,e) = - | E B n (t)cos ne f(t)dt , (2.6) 

a n 0 
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( 2 . 8 ) 


o, 9e (a,e) 



t C n (t)cos ne f(t) dt, 
n=0 


a lre( b ’ 0 ) = ir(ic+i) j ^ D n ^) sin n ® -W dt , 

f 9 • 

a l ~ ~ J ^ E n (t)cos ne f ( t) dt » (2.9) 

e n=cl 

From (2.1 )-(2.3) and (2.5)-(2.9) the Frouler coefficients may be Obtained 
as follows: 

B„(t) - C 0 (t) - - | , E 0 (t) - . | (|) , (2.10) 

A,(t) = B,(t) = 1 (f) 2 . C,(t) - - | (f) 2 , 0,(1) « -E, (t)“g£3(p) 2 -2] , 


c„(t) 

E„(t) 


( 2 . 11 ) 

J ["(I)" , -(n*2)(|) n ’) , 8„(t) = | [M)(|) n+1 *(n.E)(|) n ' 1 

a C-( n+ 2)(f) +(ti-2)(|) ], D n (t) * j- t(n+2)(g-) n '-ntp)" ' 

1 t n+1 + n-1 

FHn+2)(f) +(n+2)(i) 3, (rt = 2,3,,..) * (2.12) 


) , 


]. 


b. Mich.eJLl , & Solution 

If the loads satisfy the symmetry condition about the polar axis ox 
shown in Fig. 1, and the resultants of the tractions on the inner circu- 
lar boundary r=a are zero, then the Michel 1 * s Solution for the hollow cylinder gives 
the following stress components [4], [1] 

, f 9 b 2c. 

°2rr^ r ’ 0 ) s j 2c Q +(- -pr + 2d 1 n)cos e - 

e 


"n»2 fa n rn " 2+b n (h " 2 )rn+c n r " n " 2+d n ( n+2 > r~ n ]cosne }dt , 



a 2re^ r,e ) = j fr + 2d 1 r)s1ne 

e 


CT 200 (r * e > 

+ I [a r n “ 2 +b (n+2)r n +c_r“ n "^+d n (n-2)r" n ]cosne>dt, (2.15) 
^^2 n n n n 

where a„, b . c . and d„ (n=0,l,...) are functions of t and the constants 
n n n n 

,9 

f(t)a 

4 

e 

replace the regular coefficients in the Mitchell's solution. 

3. Boundary Conditions 

Consider now a hollow elastic cylinder having radii a and b and 

* I 

containing a radial crack along 9=0, e<r<g as shown in Figure 1. Assume 
that the inner boundary of the cylinder is reinforced by a thin mem- 
brane of thickness h and the self-equilibrating normal tractions on the 
crack surfaces are the only nonvanishing external loads. 

Let the stress state in the cracked hollow cylinder be the sum of 
the two fundamental solutions given by (2.1)-(2.3) and (2.13)-(2.15) , 
that is, let 

°ij(^»0) = * 0 » J = F»0) . (3.1) 

Then the boundary conditions of the problem may be expressed as 

a-, rr (b,e) + a 2rT (b,e) * 0, <^ r0 (M) + ^re^* 0 ^ “ °* (°l e<2ir ) 


n (t)dt , | f(t)b n (t)dt , ... 
e 


+ z [a r n “ 2 +b n n r n -c n r" n “ 2 -d„n r* n ]sin ne)dt , (2.14) 

n=s 2 n n n n 


b A 2c. 

+ 2c q + (-^3- + 6d^r)cose 
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(3.2) 



<> lrr (a>e)+o 2 rr (a,e)=p r (e) , a lr 0 (a,e)+cr 2 re (a,e)*p e (e), (O<0<2 tt) 

(3,3) 


and 


0 lee( r ’+°) + a 2ee^ r,?0 ^ = q ( r ) » ( e<r< 9) 


(3.4) 


where p r and p Q are the normal and shear tractions at the cylinder- 
reinforcement interface (Figure lb). From the equilibrium of the mem- 
brane it may easily be shown that 


P r (e) = <* rr (a»e) = xa 00 (a,e) , (O<0<2 tt) 

P fl (0) = CT r0^ a ’ 0 ^ = ' x a flfl( a ’ e ) * (O<0<2 tt) 


30 00 ' 


(3.5) 

(3.6) 


where 


hE 


— jp for plane stress , 

E m h(l-v2) 

iEn-vJ)+ht v(l+v) for p,ane stra1n - 

m m 


x = < 


(3,7) 


In (3.7) E and v are the elastic constants of the cylinder and E and v 

mm 

of the reinforcing membrane. 

By defining now 

a = . Jj. f(t) 5 _ P_f(t) ■ c 

n tt(k+ 1 ) a n 5 n ^(k+1 ) ’ 


c n " tt(k+1 ) Y n * d n = tt(k+T) V ( n=0 ^ •) 


(3.8) 


and by substituting the stresses into the boundary conditions (3.2), 
(3.5) and (3.6) we obtain 


ft 


b? + 2y o = E o ’ (1+x) a? + 2 ( ] ‘ x )y 0 = XC Q - B , 


(3.9) 
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(3,10) 


^ Y 1 Y 1 
" -p- + 26 jb = -E ] , 2(l+x) ^-+ (6x-2)d 1 a = B-j-XCj 

2Yi Y, 

- ”jjr+ 26-jb = D-j , 2(l+x) pr+ (6x-2)6 1 a * -A-j-XC-j , (3.11) 

“n‘ >n " 2+8 n (n ' 2)bn % b ’ <n+2)+s n (n+2)b ’ 1 ’ " E n (3-12) 

“n bn ' 2+s n n »"-T n b ' (n+2) -« n " k" n = D n (3.13) 

a n a n ' 2 (l+x)+6n an [(n- 2 )+x( n + 2 )]+Y n a"( n+2 *(l+x)+«na" n [(n+ 2 )+A(n-2)] 

" V XC n ( 3 -' 4 ) 

a n a n ‘ 2 0-nx)+6 n a n [ n -x n ( n +2)]- Yn a' (n+2) (l+ n x)-s n a‘ n [n+xn(n-2)] 


= A+X_C 
n n n 


(3.15) 


for n-2,3,... Note that A.., are known functions and are given by 

(2. 10)— (2.12) . Equations (3.10) and (3.11) imply that the constants A-j , 
B-j, D-j and E-j must satisfy 

V B 1 =0 > VV 0 • (3.16) 

It can be shown that equations (3.16) indeed follow from the fact that 
in the "crack solution" c^. the self-equilibrating crack surface trac- 
tions are the only non-zero external loads and consequently the stress 
vector acting on r = a has a zero resultant [1]. 

With (3.16), the solution of (3.9)-(3.11) is found to be 

8 o ■ (t - T> • <3 -’ 7) 

Y o * P'(l-x)-b2{1+x) C-0 + x)t + ] , (3.18) 
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v = a 4 (l-3x) r 3t 2 . . 2 . b^n 

Y 1 l-3X>-b^(T+A) *• 2~ + b2 + 2t^ ’ 


(3.19) 


6 1 = "a 4 (l-3x)-b^(l+; 


+ b 2 (l+x) + ^ fc 3X > ]. (3.20) 


Similarly, the solution of the infinite system (3.12)-(3.15) can be 
obtained as follows: 


a n = - e n (n-l)b2 - « b* 2 ^-D , 


6 

n a. 


(3.21) 

(3.22) 


Y n = m C-a 2(n+1) f (n+,) (l-2nx-A)-J n a2(l+n-n)i-A)+e n a 2(n+1) (l-2nx-A)] , 

A 6 (3.23) 

6 n = ^ * (3.24) 

where the determinants a q , a $ and a^ are given by 

a q = n 2 (l+2A+A 2 )a -2 b 2 -2(n 2 -l)(l-v 2 )-[l+2nx+(2n-l )A 2 ]a“ 2n b 2n 

- [l“2nx-(2n+l)x 2 ]a 2n b' 2n + n 2 (l -2x-3x 2 )a 2 b" 2 , (3.25) 

A g = {(1 -n 2 ) (1-X 2 )-[1 -2nx-(2n+l )x 2 ]a 2n b 2n +n 2 (l -2x-3X 2 )a 2 b" 2 }t"^ n+1 ) 

- + U2+n-n 2 )(l-x 2 )b- 2 -(2 + n-n 2 )(l + 2x + X 2 )a- 2 }t“ (n - 1 ) + {(l + n)[l + 2nx-(l-2n)x 2 ]a" 2n 


-0+n)(l-x 2 )b" 2n >t n ‘ 1 +{(2+n)(l+2x+x 2 )a" 2 b" 2n 
“(2+n)[l+2nx-(l-2n)x 2 ]a" 2n b" 2 >t n+1 , 


(3.26) 
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a 6 = {-(l-n)[l- 2 nx-(l+ 2 n)x 2 ]a 2n +(l-n)(l-x 2 )b 2 V"^ n+1 )-{( 2 -n)(l+ 2 xn 2 )a" 2 b 2n 
-(2-n)[l-2nX-(l+2n)x 2 ]a 2n b‘ 2 }t" (n '^+{n 2 (l+2x+x 2 )a _2 b 2 +(l-n 2 )(l-x 2 ) 
-D-2nx-(l+2n)x 2 ]a 2n b" 2n }t n “ 1 +{(2-n-n 2 )(l+2x+x 2 )a" 2 -(2-n-n 2 )(l-x 2 )b“ 2 }t n+1 . 

(3.27) 


4. The Integral Equation 

Having determined a n , B n » v n » and S n> from (3.8) and from the expres- 
sions for stresses it is seen that the solution of the cracked cylinder prob- 
lem is completely known in terms of the crack surface displacement derivative 
f(r). This function is as yet unknown and may be determined by using the 
remaining boundary condition (3.4). Thus, by substituting the stresses cr-j^ 
anc * CT 2ee as ex P resse d in terms of f(r) into (3.4), after some manipulations 
and separating the singular part of the kernel we obtain the following 
integral equation: 

g g 

t=r~ dt + / k(r,t)f(t)dt = - ^-^(r), e<r<g , (4.1) 

e e 

where the bounded kernel k(r,t) is given by 



+ n=2 Can,-n " 2 + 6 n (n+2)r " + \ r ' < ' n * Z) + <! n (n-2)r _n ]} (4.2) 

From the definition of f(r) as given by (2.4) it is seen that the 
following singlevaluedness condition must be satisfied 

f(r)dr = 0 . (4.3) 

e 
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The singular integral equation (4.1) may easily be solved numerically 
by using the procedure outlined, for example, in [7]. Noting that 
the solution of (4.1) is of the form [8] 

f(r) = F(r)//(r-e)(g-r) , ( 4 , 4 ) 

where F(r) is a bounded function, after determining f(r) the stress 
intensity factors at the crack tips may be defined and evaluated as 
follows: 

k(e) = lim /2(e-r) o 00 (r,O) = lim ^ /2(r-e) f(r) , (4.5) 

k(g) ~ l lm >^2( r-g) a flfl (r,0) = -lim /2( g-r) f(r) . (4.6) 

r+g 00 r->g 1 K 


5 • Special Cases 

The first special case to be considered is an infinite plane which 
contains a reinforced circular hole and a finite radial crack. In this 
case by letting b-*» the Fredholm kernel k(r,t) in (4.1) may be simplified 
as 


lim k(r,t) = k,(r,t) = - \ {- 
b-*» 1 2 



z 

n=2 


[Y n r " (n+2)+6 n (n ~ 2 ) r " n ^ 


> 



* Y 1 = 



1 

t* 


(5.1) 

(5.2) 


Y 


n 


a 2 (n + l)( n 2 - 3 An 2 ) } + a 2 n (n 2 -n- 2 -An z +An+ 2 A) 1 

2 An + 1 -A * t n+l 2 An+l-A * n- 1 . * 


= a 2 n (n-An-l+A) 
n 2 An+l-A 


1 



a 2 ^ n " 1 )( n - 2 +An- 2 A) 1 

2 An+l-A * t n-l 


(5.3) 

(5.4) 
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Note that the new kernel given by (5.1) can be evaluated partly in closed 
form by using the following infinite series: 


Ax + (A+B)x 2 + (A+2B)x 3 + ... = - 

where A and B are arbitrary constants. 

After some obvious transformations it can be shown that 

k f r t ) - - 1 t is. + -1 + 1 r 15a6 3a4 + a4 + 3a \ 1 

k l (r,t) 2 { F + 7^ 4 ( r 3 t 5r " ~F~ rt 2 “F' rtt 


(5.5) 


1 ,3a 6 


1 ,3a 4 A 9a 2 3a 2 

~f 


J_ a- a^ a^ a* I ^a 1 * . 

2 [ rW " F ‘ rF " r ' (rt-a 2 ) 2 " 4 T*~ ~ rt? 

2.n-l . a 4 a 2 a 2 1 . 00 n2 


+ 23 U - 1 ,aV , a 4 a 2 a 2 . 1, - n 2 

E 2xn+i-x ( FF - ( FF- ' F ‘ Ft + P i 0 2IFF 


n=2 

n-1 


n=2 

n-1 


,a 2 i 1,3a 4 . a 2 . 5a 2 9x - n ,a 2 x n_1 . 1,3a 4 

rt' - 2 ( FF + F + H 2 " PJ, 2xFTF rt } + 4 ( FF 


. 9a 2 3a 2 9* • 

+ IF" " Ft* ' P n ? 2 


* (§1.) 

2xn+l-x W 


n=2 
2 n-1 


(5.6) 


Now observing that 

oc 1 a2 

1 (Sr) 


n-1 


rt 


and 


n=2 2xn+l-x 'rt ; 2xa 2 ^ n ^rt^ 


l0 9 <T^> * f T • ■ 1 i x<1 > 


1 fSi," _ JUzil 


izii ; i ,«£> 

2x n _2 n(2nx+l-x) 'rt' 


n-1 


(5.7) 


(5.8) 


the kernel k-|(r,t) may be decomposed as follows: 
k^r.t) = k ls (r,t) + k lws (r,t) + k ]b (r,t) 


(5.9) 


where 
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where 



■ - i 



1 

rt-a 7 


1 r 3a6 .. _ a4 _ a4 af_\ a 2 

4 'r 3 t 2 r 7 ~ rt 2 " r (rt-a 2 ) 2 * 


rt /3a 4 


i, 4 * \ - rt /Od’ , 9a 2 ja fc 

k iws (r,t) " ■ TePx ( Ft 2 ' + 73" Tt 7 


3a 2 ^ 23 


)log(rt-a 2 ) 


(5.10) 

(5.11) 


1 / (r ti - 1 r 3a4 + 9a2 3a 2 , 23\ r rt 1on 

k lb (r ’ t} ' 8 + W ‘ Ft 2- + T )C 2Px 109 rt 


X 

2x 


n=2 


(i-*)({£) n '’ 

2xn(2xn+l-x) 

■2ti 


1 r p o . ‘ Y 1 /a 4 a 2 a 2 , 1 x “ 

- 2 { " fJ- + -pr " W ■ W ft 7 + F> „ E , 

n=2 


n /aM-l 

1 / ja- + 11 + 211 _ 1 \ y _}rV 

2 ^r 3 t 2 r 7 rt 7 r ' _ 2 2xri+l-x 


1 / 3a 4 . a 2 . 5a 2 9 


*<Ft> 


2 n-1 


, 1 /3a 4 . 9a 2 3a 2 9x » 

+ 4 'Ft* + -W - ■ P n f 2 2xn+l-x 


2xri+l -x 


(5.12) 


The subscripts in k-j $ , k-j ws and k^ refer to the fact that as r and t both 
approach e = a (in the case of the crack touching the interface), these 
kernels become "singular", "weakly singular" and "bounded", respectively. 

The second special case is that of an infinite plane containing 
a radial crack and a circular hole without reinforcement. In this case 
by letting X = 0 (see (3.7)) it can be shown that the kernel k(r,t) in 
(4.1) may be evaluated in closed form as follows: 


lim k(r,t) = k 2 (r,t) = - ffit-affi 2 + 

x->0 


t(t 2 -a 2 ) 

(rt-a 2 ) 2 


+ 


t . a 2 -t 2 _ a 2 1 
rt-a 2 r 2 t r 7 " r ‘ 


(5.13) 
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The third special case refers to the infinite plane with a rein- 
forced circular hole and a radial crack in which the radius a is 
infinitely large. This is the problem of a cracked half plane rein- 
forced by a membrane along its boundary (Figure 2). In this case if we 
defi ne 

r = a + y, t = a + 5»e = a + e 0 , g = a + g 0> (5.14) 

and assume that y and 5 remains finite as a tends to infinity, the 
kernel k(r,t) in (4.1) becomes 

«»-.*) | = k](a+y,a+c) « k 3 (y, ? ), (e o <(y, ? )<g o ) . (5.15) 

a-*» 


Again, by evaluating the asymptotically dominant terms in the infinite 
series in closed form, after somewhat lengthy manipulations we obtain 


k 3 (y,t) = - c ] y - *hf|+y) “ 2^h%iy7 + iFfc+y ) 2 


+ (4 + 2^ + ^>- ^ j 2 fW^rr 


n -1 


(5.16) 


where 


♦ - 


E m /E for plane stress, 
E m 0-v 2 )/E(l-v m 2 ) for plane strain . 


(5.17) 


Of course, the integral equation for the half plane problem may also 
be derived directly by using the method of Fourier transforms. By follow 
ing this technique and by omitting the details, after evaluating the 
asymptotically dominant infinite integral si the kernel k 3 (y,c) may be 
obtained as follows: 



(5,18) 


- 13 - 



Finally, for the half plane without reinforcement it can be shown 
that for a-*» (5.13) would reduce to the following known kernel [7] 

k 4 (y,?) ' ' + Tc+yF ' Ts+yF • < 5 - 19 > 


6. The Case of the Crack Touching The Interface (e=a) 

In the infinite plane containing a reinforced circular hole considered 
in the previous section, let us assume that the inner crack tip touches 
the interface, (i.e., e=a). In this case, from (5.9)-(5,12) the kernel 
k(r,t) of the integral equation (4.1) may be shown to reduce to 


lim k(r,t) = k^r.t) = k ]s (r,t) + k lws (r,t) + k lb (r,t) 


( 6 . 1 ) 


e->a 


/ r t > - 1 (3a*» 6a2 ^ a 2 (r+a) 2 (r-a) 2 1 

8 ( F F' 5 ' 7 a£ 4r 5 ' * “ 

r " r 

k lws (r>t) = ■Wx ( f^ + ^-^ + 7 ) log(rt-a 2 ) 


tt- ( 6 - 2 ) 
(t - ~) 2 


(6.3) 


ki b (^t) = k lb (r,t) + + 5) + jk ($£ + M). 


‘lb' 


8t ' r 


8t^ v r ; 


(6.4) 


Note that as both r and t approach the interface r = a, the kernel 
k 1$ become unbounded and hence influence the singular behavior of the 
solution f(t). This behavior may be investigated by letting 

f(t) = (g-t)" a (t-a)“ 6 , (0<Re(a,e)<l ) (6.5) 

and by using the function theoretic method described in [7]. Omitting 
the details, the characteristic equations giving a and e may thus be 
obtained as 

cos to - 0 , cot ire - s -- .^ g = 0 , (6.6) 
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giving a - 1/2 and $ = 0 as the acceptable roots. The root e = 0 indicate 
that at the crack tip r = a the crack surface derivative is finite. 

Since the bending stiffness of the reinforcement is zero, this is the 
physically expected result. 

In the limiting case of a-*° (Figure 2), the kernel given by (5.18) 
which is obtained from the method of Frouier transforms is still valid. 

In this case the integral in (5.18) may be evaluated as 

r e -(s + y)<* -0 00 A n 

J 2 *h« + l d “ ■ - Sfh [e o +1 °9(|^) + ! HTprr] • ( 6 - 7 ) 

0 1 

where 0 = -(?+y)/2<frh and E Q = 0.5772157 is the Euler's constant. Hence 
the infinite integral in (5.18) has the following order: 

f e -« (c+y) 

J 2 <j>ha+l da ~ o nog(c+y)] (6.8) 

0 


and the kernel of the integral equation (4.1) may be arranged as follows: 

k(r,t) = k 3 (y.s) = C 3s (y.c) + ^ 3ws (y.c> + k 3b (y.c) , (6.9) 


W*’ 4 ’ ' - th ' 
< 

k 3ws ( y-'> ■ 4 

4 



da , 


( 6 . 10 ) 

( 6 . 11 ) 



It is seen that as 5 and y both approach zero k 3s becomes unbounded, ic^ 
exhibits only a weak singularity and k 3fa remains unbounded. 
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Thus, if we express f(t) again as in (6.5), the function theoretic 
method gives the characteristic equations for a and $ as 

cos uo - 0 , (cosTr6-l)/sinir6 = 0 (6.13) 

which are the same as (6.6) and give a = 1/2 and 6 = 0. This result 
could also be obtained from the circular hole problem by substituting 

r “ a+ y» * 58 a+ 5» an£ * by letting a-*» in (6.2). After some simple algebra 
one may then observe that the principal part of the kernel k ls is 
-V(c+y) which is identical to (6.10). S 


7. Solution and Results 

Excepting the problem of unreinforced edge crack, in all cases con- 
sidered the displacement u Q at the crack tips is zero and consequently 
the integral equation (4.1) and all its special cases must be solved under 
the single-valuedness condition (4.3). For the numerical solution follow- 
ing normalized quantities are introduced: 


t = ^ T + Sje , r = ^ p + 3je , 


(7.1) 


k(p»x) = -9^- k(r,t) , q(p) = q(r) , 
f(x) = f(t) = F(t)w(t) , (— 1<T<1) 


(7.2) 

(7.3) 


In the embedded crack problem (i.e., for e>a) w(t) = (1- t 2 )"^ 2 and 
a Gauss-Chebyshev Integration technique is used to determine F(t) [7]. 

The stress intensity factors are then evaluated from (see (4.5) and (4.6)) 


k(e) - /(g-ej/2 F(-l) , k(g) = -/(g-e)/2 F(l). (7.4) 


For e=a the function f(r) is bounded at r=a, (i.e. in (6.5) 6=0) and the 
stress intensity factor at r=g is obtained from 
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k(g) = -/Fa F(l) . 


t7.5) 


The calculated stress intensity factors are given in Tables 1-6. 
In the case of an embedded crack (i.e., for e>a or e Q >0) the stress 
intensity factors are normalized as follows: 


— k < e > -Me) 

■,9 0 /(g-e|/2 


= k (q) 




(7.6) 


Table 1 gives an example showing the effect of the stiffness parameter 
A defined by (3.7) on the stress intensity factors. Similar results 
for a reinforced half plane are shown in Table 2. The stiffness con- 
stant <f> shown in Table 2 is defined by (5.17). 

The results for somewhat more realistic values of the stiffness 

parameter A and the radius ratio a/b which may simulate a pressure vessel 

with internal cladding are shown in Table 3-6. The normalized stress 

intensity factors are calculated for two types of external loads. In 

the first it is assumed that a uniform crack surface pressure a (r,0) = 

/ \ • 

( Hn = -q 0 is the only nonzero applied load. In the second the cylinder 
is assumed to be under internal pressure a rr (a-h,e) = -p Q . 

Generally, for the embedded crack the influence of the membrane 
reinforcement or the cladding on the stress intensity factors does not 
seem to be very significant. Since the properties of the clad material 
cannot be very different than that of the cylinder, this result is not 
surprising. However, in the case of the crack terminating at the inter- 
face (i.e., for e=a, Figure 1) the stress intensity factors k(g) cal- 
culated for Aft) and A=0 are very different. For A=0 there is no reinforce- 
ment and the crack is an edge crack with a non-zero opening at r=a 
(i.e., u 0 (a,+O)-u 0 (a,-O)>O). On the other hand, for A^O since u 0 (a,O)=O 
there is a great deal of constraint tending to reduce the stress intensity 
factor. One may note that because of the membrane assumption made for 
the reinforcement at r=a the derivative of the crack surface displacement 
is finite giving a finite angle between the upper and lower crack surfaces. 
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Table 1. Stress intensity factors in a reinforced hollow circular 
cylinder containing a radial crack (Figure 1). External 
load: uniform crack surface pressure q(r) = -q , a = 0.5b, 
g-e - 0.5(b-a), k n (e) = k(e)/q Q ,/(g-e)/2, k n (g) = k(g)/q 0 /(g-e)/2. 


e-a 

X 

= 0 

X 

= 0.2 

X = 

0.5 

b-a 

k n (e) 

kn(g) 

k n< e ) 


M e > 

Ms) 

0.15 

1.1968 

1.1754 

1.1085 

1.1255 

1.0549 

1.0988 

0.20 

1.1517 

1.1769 

1.1123 

1.1436 

1.0740 

1.1223 

0.25 

1.1331 

1.1951 

1.1215 

1.1706 

1.0926 

1.1523 

0.30 

1.1727 

1.2303 

1.1379 

1.2100 

1.1148 

1.1930 

0.35 

1.1893 

1.2912 

1.1637 

1.2702 

1.1439 

1.2526 

0.40 

1 .2263 

1.4033 

1 . 2042 

1.3764 

1.1854 

1.3559 
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Table 2. Stress intensity factors in a reinforced half plane containing 
a crack perpendicular to the boundary (Figured). External 
load: uniform crack surface pressure q(y) = -q rt , = h<t/q 4 

W = k (9 0 )/V ( VV /2 • k n (e o> = k ( e o)^°#y7?. ' 


e 0 /9o 

X o 

= 0 

X o = 

0.0375 

X o = 

0.075 

X o = 

0.15 

x o = 

0.3 

k n< e o> 

k n<9 0 ) 

W 

W 

k n (e o> 

k n(9o> 

k n< e o> 

k n<9n) 

k n ( e 0 ) 

k n(9 0 ) 

0.1 

1.42667 

1.15500 

1.41438 

1.15234 

1 .40300 

1.14983 

1 .38222 

1.14510 

1.34633 

n u 
1.13661 

0.2 

1 .20354 

1.09670 

i 

1.19970 

1.09543 

1.19604 

1.09422 

1.18921 

1.09190 

11.17698 

1.08765 

0.3 

1.11075 

1.06242 

1.10909 

1.06171 

1.10750 

1.06103 

1.10449 

1.05972 

1.09902 

1.05730 

0.4 

1.06224 

1.03990 

1.06145 

1.03949 

1.06068 

1 .03909 

1.05923 

1.03833 

1.05655 

1.03691, 


Table 3. Stress Intensity factors in a reinforced hollow circular cylinder 
containing a radial crack (Figure 1). External load: uniform 
crack s urface pr essure q(r) = -q o , a/b = 0.8, e/a = 1.05, k (g) = 
k(g)/q 0 vT9-e)/2 , k n (e) = k(e)/q /fg-e)/4 . " 


h-a 

O 

11 

c< 

A = 0.001 

X = 

0.002 

X = 

0.004 


k n (e) 

k „(9) 

k n (e) 

k n(g> 

k n (e) 

k n (g) 

k „(e) 

k n (g) 

0.05 

1.0044 

1.0041 

1.0043 

1.0039 

1.0042 

1.0038 

1.0039 

1.0036 

0.10 

1.0157 

1.0133 

1.0153 

1.0129 

1.0148 

1.0126 

1.0140 

1.0119 

0.20 

1.0511 

1.0382 

1.0497 

1.0373 

1.0484 

1.0364 

1.0459 

1 .0349 

0.30 

1.0962 

1.0668 

1.0938 

1.0656 

1.0915 

1.0644 

1.0871 

1 .0621 

0.40 

1.1466 

1.1023 

1.1432 

1.1011 

1.1400 

1.0998 

1.1339 

1 .0974 

0.50 

1.2021 

1.1596 

1.1979 

1.1584 

1.1940 

1.1573 

1.1866 

1.1552 
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Table 4. Stress intensity factors in a reinforced pressurized cylinder 
containing a radial crack. External load: internal pressure 

a rr = -P 0 » a / b = °* 8 » e / a = k n ( e ) = k (e)/P 0 /(g-e)/2 , 

k n (g) = k(g)/p 0 Ag-e)/2 . 



X « 



X = 0.004 

b-a 

k n (e) 

k„(g) 

k„(e) 

k„(g) 

k„(e) 

k„(g) 

0.05 

4.3014 

4.2703 

4.3003 

4.2693 

4.2992 

4.2684 

0.10 

4.3346 

4.2662 

4.3307 

4.2633 

4.3273 

4.2606 

0.20 

4.4534 

4.2872 

4.4420 

4.2800 

4.4316 

4.2734 

0.30 

4.6101 

4.3256 

4.5903 

4.3156 

4.5723 

4.3064 

0.40 

4.7854 

4.3935 

4.7581 

4.3828 

4.7332 

4.3730 

0.50 

4.9779 

4.5473 

4.9451 

4.5382 

4.9150 

4.5297 


Table 5. The stress intensity factor k n (g) = k(g)/q Ql /g^a" in a reinforced 
cylinder containing a radial crack which terminates at the 
interface (e=a, Figure 1). External load: crack surface pressure 
q(r) = -q Q » a/b = 0.8. 


S=§L 

b-a 

II 

O 

X = 0.002 

X = 0.004 

0.25 

1.3380 

0.8894 

0.8701 

0.30 

1.4477 

0.9256 

0.9023 

0.40 

1.7159 

0.9886 

0.9580 

0.50 

2.0489 

1.0375 

1.0013 

0.60 

2.4352 

1.0735 

1.0345 

0.70 

2.8271 

1.1065 

1.0681 
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Table 6. Stress intensity factor k n (g) = k(g)/p o ^iT in a pressurized 
reinforced cylinder containing a radial crack which terminates 
at the interface (e=a, Figure 1). External load: o rr (a,0) = 
-p Q , a/b = 0.8. (For the edge crack case, that is, for x = 0 
the effect of the pressure a 00 (r,O) = -p Q on the crack surfaces 
is not included). (k n *(g) * k n (g)/o 00 (a,O)^ 


b-a 

X = 0 

k n (g) 

X » 0 

k„*(g) 

X = 0.002 
k n (9 > 

X . 0.004 

k„(g) 

0.25 

5.8477 

1.2836 

3.8533 

3.7684 

0.30 

6.2838 

1.3794 

3.9755 

3.8734 

0.40 

7.3578 

1.6151 

4.1768 

4.0443 

0.50 

8.6933 

1.9083 

4.3154 

4.1603 

0.60 

10.2352 

2.2467 

4.3977 

4.2320 
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Figure 1. Geometry of a reinforced thick-walled cylinder containing 
a radial crack. 
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Figure 2. 




Notation for an elastic half plane which is reinforced by a 
membrane and which contains a finite crack. 
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